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Abstract 

We consider a dynamic vehicle routing problem in wireless networks where messages arriving randomly in time and space 
are collected by a mobile receiver (vehicle or a collector). The collector is responsible for receiving these messages via wireless 
communication by dynamically adjusting its position in the network. Our goal is to utilize a combination of wireless transmission 
and controlled mobility to improve the delay performance in such networks. We show that the necessary and sufficient condition 
for the stability of such a system (in the bounded average number of messages sense) is given by p < 1 where p is the average 
system load. We derive fundamental lower bounds for the delay in the system and develop policies that are stable for all loads 
p < 1 and that have asymptotically optimal delay scaling. Furthermore, we extend our analysis to the case of multiple collectors 
in the network. We show that the combination of mobility and wireless transmission results in a delay scaling of Q(j^) with 
the system load p in contrast to the €>( rrj^p ) delay scaling in the corresponding system where the collector visits each message 
location. 

I. Introduction 

There has been a significant amount of interest in performance analysis of mobility assisted wireless networks in the last 
decade (e.g., Q4), (20|, EH, J26), E3, (29|, ED, E3). Typically, throughput and delay performance of networks were 
analyzed where nodes moving according to a random mobility model were utilized for relaying data (e.g., |fl3ll , lfl4l . |23ll ). 
More recently, networks deploying nodes with controlled mobility have been considered focusing primarily on route design and 
■ ignoring the communication aspect of the problem (e.g., ||20l , ||29l , |30l , l33l ). In this paper we explore the use of controlled 
mobility and wireless transmission in order to improve the delay performance of wireless networks. We consider a dynamic 
vehicle routing problem where a vehicle (collector) uses a combination of physical movement and wireless reception to receive 
randomly arriving data messages. 

Our model consists of a collector that is responsible for gathering messages that arrive randomly in time at uniformly 
distributed geographical locations. The messages are transmitted when the collector is within their communication distance 
and depart the system upon successful transmission. The collector adjusts its position in order to successfully receive these 
messages in the least amount of time as shown in Fig. Q] This setup is particularly applicable to networks deployed in a large 
area so that a mobile element is necessary to provide connectivity between spatially separated entities in the network. For 
instance, this model is applicable to a sensor network where a mobile base station collects data from a large number of sensors 




Fig. 1. The collector adjusts its position in order to receive randomly arriving messages via wireless communication. The circles with radius r* represent 
the communication range and the dashed lines represent the collector's path. 



deployed at random locations inside the network J9), l32l , l33ll . Another application is utilizing Unmanned Aerial Vehicles 
(UAVs) as data harvesting devices or as communication relays on a battlefield environment l25ll . This model also applies to 
networks in which data rate is relatively low so that data transmission time is comparable to the collector's travel time, for 
instance in underwater sensor networks [fl~), l28l . 

Vehicle Routing Problems (VRPs) have been extensively studied in the past (e.g., flU, |6), Q, (8), 0J], ED, ED, ED, 
||3T1 ). The common example of a VRP is the Euclidean Traveling Salesman Problem (TSP) in which a single server is to 
visit each member of a fixed set of locations on the plane such that the total travel cost is minimized. Several extensions of 
TSP have been considered in the past such as stochastic demand arrivals and the use of multiple servers @, 0, ifTTl . In 
particular, in the TSP with neighborhoods (TSPN) problem the vehicle is to visit a neighborhood of each demand location J4), 
ll2TI . which can model a mobile collector receiving messages from a communication distance. A more detailed review of the 
literature in this field can be found in Q, ETl and ll3TI . 

Of particular relevance to us among the VRPs is the Dynamic Traveling Repairman Problem (DTRP) due to Bertsimas and 
van Ryzin (6), Q, 0. DTRP is a stochastic and dynamic VRP in which a vehicle is to serve demands that arrive randomly in 
time and space. Fundamental lower bounds on delay were established and several vehicle routing policies were analyzed for 
DTRP for a single server in j6), for multiple servers in J7) and for general demand and interarrival time distributions in (8). 
Later, [29], ll30l generalized this model to analyze Dynamic Pickup and Delivery Problem (DPDP) where fundamental bounds 
on delay were established. We apply this model to wireless networks where the demands are data messages to be transmitted to 
a collector which is capable of wireless communicatiorfj In our system the problem has considerably different characteristics 
since in this case the collector does not have to visit message locations but rather can receive the messages from a distance 
using wireless communication. The objective in our system is to effectively utilize this combination of wireless transmission 
and controlled mobility in order to minimize the time average message waiting time. 

In a closely related problem where multiple mobile nodes with controlled mobility and communication capability relay the 
messages of static nodes, ll26l derived a lower bound on node travel times. Message sources and destinations are modeled 
as static nodes in 11261 and these nodes have saturated arrivals hence queuing aspects were not considered. In an independent 
work, 1171 considered utilizing mobile wireless servers as data relays on periodic routes and applied various delay relations 
from Polling models to this setup. A mobile server harvesting data from spatial queues in a wireless network was considered 
in 11251 where the stability region of the system was characterized using a fluid model approximation. In ifTOl we analyzed a 
model similar to the current paper but for which the arriving messages were transmitted to the collector using a random access 
scheme, creating interference among neighboring transmissions. In this paper, the message transmissions are scheduled, i.e., 
there is only one transmission in the system at a given time, and the collector decides on the message to be transmitted next. 
The two systems have considerably different characteristics as will be explained in the following sections. 

Another related body of literature lies in the area of utilizing mobile elements that can control their mobility to collect 
sensor data in Delay Tolerant Networks (DTN) (e.g., 0, ED, ED, ED, E2, G3)- Route selection (e.g., ED, El, ED), 
scheduling or dynamic mobility control (e.g., (9), 11281 . ll32l ) algorithms were proposed to maximize network lifetime, to 
provide connectivity or to minimize delay. These works focus primarily on mobility and usually consider particular policies 
for the mobile element. To the best of our knowledge, this is the first attempt to develop fundamental bounds on delay in a 
system where a collector is to gather data messages randomly arriving in time and space using wireless communication and 

'in 1 6 1, (7J, or Qi] the collector needs to be at the message location in order to be able to serve it, therefore, we will refer to the DTRP model as the system 
without wireless transmission. 
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controlled mobility. 

The main contributions of this work are the following. We show that p < 1 is the necessary and sufficient condition for 
the stability of the system where p is the system load. We derive fundamental lower bounds on delay and develop algorithms 
that are asymptotically within a constant factor of the lower bounds. We show that the combination of mobility and wireless 
transmission results in a delay scaling of 9(1/(1 — p)) in contrast to the 0(1/(1 — p) 2 ) delay scaling in the system where the 
collector visits each message location analyzed in [6), JT). Finally, we extend our analysis to the case of multiple collectors in 
the network. 

This paper is organized as follows. In Section [TT] we describe the system model and in Section [Til] we characterize the 
necessary and sufficient conditions for the stability of the system. We derive fundamental lower bounds on delay in Section 
IIVI and in Section [VI we provide upper bounds on delay together with numerical results. We extend the analysis for a system 
with multiple collectors in Section [Vll 

II. Model 

Consider a square region 7Z of area A and messages arriving into 1Z according to a Poisson process (in time) of intensity 
A. Upon arrival the messages are distributed independently and uniformly in 1Z and they are to be gathered by a collector via 
wireless reception. An arriving message is transmitted to the collector when the collector comes within the reception distance 
of the message location and grants access for the message's transmission. Therefore, there is no interference power from the 
neighboring nodes during message receptions. We assume that the transmit power Pt is constant and that the transmissions 
are subject to distance attenuation. In such a system, the received power of a transmission from node i, located at distance 
from the collector, is given by Pr $ = PxKr~ a lfT4l . 03), where a is the power loss exponent (typically between 2 and 6), 
and K is the attenuation constant normalized to 1. 

Next we argue that the Signal to Noise Ratio (SNR) packet reception model 1141 . |[f5l is equivalent to a disk model iflOl . 
|[T5l under the above assumptions. In the SNR model, a transmission is successfully decoded at the collector if its SNR is 
above a threshold (3, i.e., if SNR; = Pr.i/Pn > P, where Pjv is the background noise power. We let r be the reception 
distance of the collector. If the location of the next message to be received is within r, the collector stops and attempts to 
receive the message. Otherwise, the collector travels towards the message location until it is within a distance r away from the 
message. A transmission at distance r to the collector is successful if r < (SNR,;//?) 1 /" where SNR C = Pt/Pn denotes the 
SNR of a transmission from unit distance. Therefore, the optimal reception distance is the maximum reliable communication 
distance r* = (SNRc//}) 1 /". Hence, essentially we have a disk model of radius r*, where a transmission can be received only 
if it is within a disk of radius r* around the collector. Under this model, transmissions are assumed to be at a constant rate 
taking a fixed amount of time denoted by s. 

The collector travels from the current message reception point to the next message reception point at a constant speed v. 
We assume that at a given time the collector knows the locations and the arrival times of the messages that arrived before this 
time. The knowledge of the service locations is a standard assumption in vehicle routing literature (4), J6), ifTTI . Ell . l29l . 
Let N(t) denote the total number of messages in the system at time t. We say that the system is stable under a policy if JJJ, 

nm 

limsupE[iV(i)] < oo, (1) 

t— ^oo 

namely, the long term expected number of messages in the system is finite. Let p = As denote the load arriving into the system 
per unit time. For stable systems, p denotes the fraction of time the collector spends receiving messages. 
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We define Tj as the time between the arrival of message i and its successful reception. Tj has three components: Wa,i, the 
waiting time due to collector's travel distance from the time message i arrives until it gets served, W s ,i, the waiting time due 
to the reception times of messages received from the time message i arrives until it gets served, and s, reception time of the 
message. The total waiting time of message i is denoted by W\ = Wd,i + W s ,i, hence W% = — s. The fraction of time the 
collector spends receiving messages is denoted by p, and for stable systems p = As. We let di be the collector travel distance 
from the collector's reception location for the message served prior to message i to collector's reception location for message 
i. The time average per-message travel distance of the collector, denoted by d, is defined by an expectation in the steady state 
given by d = limi_ J . 00 E[di}. The time average delays T, W, Wd and W s are defined similarly to have T = Wd + W s + s where 
all the limits are assumed to exist. T* is defined to be the optimal system time which is given by the policy that minimizes T. 

III. Stability 

In this section we characterize a necessary and sufficient condition for the stability of the system. 

A. Necessary Condition for Stability 

Theorem 1: A necessary condition for the stability of any policy is p < 1. Furthermore, the time average waiting time 
satisfies 

Proof: We first show that the unfinished work and the delay experienced by a message in the system stochastically 
dominates that in the equivalent system with zero travel times for the collector. 

Lemma 1: The steady state time average delay in the system is at least as big as the delay in the equivalent system in which 
travel times are considered to be zero (i.e., v — oo). 

Proof: Consider the summation of per-message reception and travel times, s and di, as the total service requirement of 
a message in each system. Since di is zero for all i in the infinite velocity system and since the reception times are constant 
equal to s for both systems, the total service requirement of each message in our system is deterministically greater than that 
of the same message in the infinite velocity system. Let D\, D2, D n and D 1 , D 2 , D n be the departure instants of the 
I s ', 2 nd and similarly the n th message in the original and the infinite velocity system respectively. Similarly let Ai,A%, A n 
be the arrival times of the 1 st , 2 nd and the n th message in both systems. We will use induction to prove that Di > D i for all 
i. Since the service requirement of each message is smaller in the infinite velocity system, we have D\ > D 1 . Assuming we 
have D n > D n , we need to show that D n+ i > D n+1 . 

A n +i < D n+ i - s, (3) 

hence the n + 1 th message is available before the time D n+ i — s. We also have 

D n < D n < D n+1 - s. 

The first inequality is due to the induction hypothesis and the second inequality is because we need at least s amount of 
time between the n th and n + 1 th transmissions. Hence the collector is available in the infinite velocity system before the 
time D n+ i — s. Combining this with d2Tb proves the induction. Now let D{t) and D (t) be the total number of departures 
by time t in our system and the infinite velocity system respectively. Similarly let N(t) and N (t) be the total number of 
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messages in the two systems at time t. Finally let A(t) be the total number of arrivals by time t in both systems. We have 
A{t) = N(t) - D(t) = N'(t) - D'(t). From the above induction we have D(t) < D' (t) and therefore 

N(t) > JV'(t). 

Since this is true at all times, we have that the time average number of customers in the system is greater than that in the 
infinite velocity system. Finally using Little's law proves the lemma. ■ 
Since the infinite velocity system behaves as an M/G/l queue (an M/G/l queue is a queue with Poisson arrivals, general i.i.d. 
service times and 1 server and an M/D/l queue has constant service times), the average waiting time in this system is given by 
the Pollaczek-Khinchin (P-K) formula for M/G/l queues |5] p. 189], i.e., As 2 / (2(1 — As)). Therefore we have (f2l). Furthermore, 
a direct consequence of this lemma is that a necessary condition for stability in the infinite velocity system is also necessary 
for our system. It is well-known that the necessary (and sufficient) condition for stability in the M/G/l systems is given by 
p < 1 (see e.g., or 02)). ■ 

B. Sufficient Condition for Stability 

Here we prove that p < 1 is a sufficient condition for stability of the system under a policy based on Euclidean TSP with 
neighborhoods (TSPN). TSPN is a generalization of TSP in which the server is to visit a neighborhood of each demand location 
via the shortest path |4|, 12TI . In our case the neighborhoods are disks of radius r* around each message location. TSPN is 
an NP-Hard problem such as TSP. Recently, lOTI proved that a Polynomial Time Approximation Scheme (PTAS) exists for 
TSPN among fat regions in the plane. A region is said to be fat if it contains a disk whose size is within a constant factor of 
the diameter of the region, e.g., a disk, and a PTAS belongs to a family of (1 + e)-approximation algorithms parameterized 
by e > 0. 

1) TSPN Policy: Assume the system is initially empty (at time to = 0). The receiver waits at the center of 7Z until the first 
message arrival, moves to serve this message and returns to the center. Let time t\ be the time at which the receiver returns to 
the center. At t\, if the system is empty, the receiver repeats the above process and we define t% similarly. If there are messages 
waiting for service at time t\, the receiver computes the TSPN tour (e.g., using the PTAS in ED ) through all the messages 
that are present in the system at time t\, receives these messages in that tour and returns to the center. We let t-z > t\ be the 
first time when the receiver returns to the center after receiving all the messages that were present in the system at t\ and 
repeat the above process. We define the epochs ti as the time the receiver returns to the center after serving all the messages 
that were present in the system at time ij_i 0. 

Let the total number of messages waiting for service at time ti, JVj = JV(ij), be the system state at time ti. Note that JVj is 
an irreducible Markov chain on countable state space N. We show the stability of the TSPN policy through the ergodicity of 
this Markov chain. 

Theorem 2: The system is stable under the TSPN policy for all loads p < 1. 
Proof: Given the system state JVj at time U, we apply the algorithm in lETl to find a TSPN tour of length L; through the 
N neighborhoods that is at most (1 + e) away from the optimal TSPN tour length L*. Note that L* can be upper bounded 
by a constant L for all JVj. This is because the collector does not have to move for messages within its communication range 
and a finite number of such disks of radius r* can cover the network region for any r* > 0. The collector then can serve the 

2 A similar policy based on TSP was discussed in 1241 for a system without communication capability similar to DTRP, where the time average delay of 
the policy was characterized for the heavy load regime. 
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messages in each disk from its center incurring a tour of constant length L (an example of such a tour is shown in Fig. |2}. We 
will use the Lyapunov-Foster criterion to show that the Markov chain described by the states Ni is positive recurrent Q. We 
use V(iVj) = sNi, the total load served during i th cycle, as the Lyapunov function (note that V(0) = 0, Sk = {x : V(x) < K} 
is a bounded set for all finite K and V(.) is a non-decreasing function). Since the arrival process is Poisson, the expected 
number of arrivals during a cycle can be upper-bounded as follows: 

E[A l+1 |A ?; ] < X(L/v + sNi). (4) 

Hence we obtain the following drift expression for the load during a cycle. 

E[sN i+1 - sNi\Ni] < pL/v - (1 - P )sN l . 
Since p < 1, there exist a 5 > such that p + 8 < 1: 

E[sA l+ i - sN^N,} < pL/v - 5sN q 

< .l {JVl6 s}, (5) 

v 

where ljjves} is equal to 1 if N £ S and zero otherwise and S = {N € N : N < K} is a bounded set with K = [A- + 1]. 
Hence the drift is negative as long as JVj is outside a bounded set. Therefore, by the standard Lyapunov-Foster criterion J2|> 
(3), the Markov chain (iVj) is positive recurrent, it has a unique stationary distribution and we can bound the steady state time 
average of Ni as |22l 

IimsupE[JV(ti)] < AL (6) 

tj-KX) U(l — p) 

Furthermore, given some t € we have 

limsupE[iV(t)] < limsupE[iV(ii) + N{t. l+1 )} 

< 2 n XL . < oo. (7) 

■ 

The delay scaling of the TSPN policy with load p is as shown in ©, the same delay scaling as in a G/G/l queue. This 
is a fundamental improvement in delay due to the communication capability as the system without wireless transmission in 
(6) has 6( ^ 1 _ 1 p p ) delay scaling. 

Note that p < 1 is a sufficient stability condition also for the system without communication capability. This case corresponds 
to r* = 0, where we utilize a (1 + e) PTAS for the optimal TSP tour through the message locations instead of the TSPN tour. 
An upper bound on the TSP tour for any A^, points arbitrarily distributed in a square of area A is given by ^/2lDVi + 1.75v / A 
030- Similar arguments as above leads to the drift condition 



E[sA i+1 - sNi\Ni] < p( Kl ^AN~+K 2 ) - (1 - p)sN ll 

for some constants k% and K2, where the drift is again negative as long as Ni is outside a bounded set S. The difference in 
this case is that the travel time per cycle scales with the number of messages N as y/Nl which can be shown to result in 
0( (il p yi ) delay scaling with the load p. 
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IV. Lower Bound On Delay 



For wireless networks with a small area and/or very good channel quality such that r* > J A/2, the collector does not need 
to move as every message will be in its reception range if it just stays at the center of the network region. In that case the 
system can be modeled as an M/D/l queue with service time s and the associated queuing delay is given by the P-K formula 
for M/G/l queues, i.e., W = As 2 /(2(1 — p)). However, when r* < \J A/2, the collector has to move in order to receive some 
of the messages. In this case the reception time s is still a constant, however, the travel time per message is now a random 
variable which is not independent over messages (for example, observing small travel times for the previous messages implies 
a dense network, and hence the future travel times per message are also expected to be small). Next we provide a lower bound 
similar to a lower bound in |6] with the added complexity of communication capability in our system. 

Theorem 3: The optimal steady state time average delay T* is lower bounded b^ 

. E\(\\U\\ - r*)+] As 2 

t *>^n — -LA + — - + s. (8) 

v{l-p) 2(1 -p) 

Proof: As outlined in Section [II] the delay of message i, Ti has three components: Ti = Wd,i + W s .i + Si. Taking 
expectations and the limit as i — > oo yields 

T = W d + W S + s. (9) 

A lower bound on Wd is found as follows: Note that Wd,i-V is the average distance the collector moves during the waiting time 
of message i. This distance is at least as large as the average distance between the location of message i and the collector's 
location at the time of message z's arrival less the reception distance 7-*. The location of an arrival is determined according to 
the uniform distribution over the network region, while the collector's location distribution is in general unknown as it depends 
on the collector's policy. We can lower bound Wd by characterizing the expected distance between a uniform arrival and the 
best a priori location in the network that minimizes the expected distance to a uniform arrival. Namely we are after the location 
v that minimizes E[||J7 — where U is a uniformly distributed random variable. The location v that solves this optimization 
is called the median of the region and in our case the median is the center of the square shaped network region. Because the 
travel distance is nonnegative, we obtain the following bound on Wd' 

Wd > SEl^H. (1 o) 

v 

Let A?" be the average number of messages received in a waiting time and let R be the average residual reception (service) 
time. Due to the PASTA property of Poisson arrivals (Poisson Arrivals See Time Averages) (see for example J5] p. 171]) a 
given arrival in steady state observes the time average steady state occupancy distribution. Therefore, the average residual time 
observed by an arrival is also R and is given by Xs 2 /2 J5J p. 188] and we have 

W s = sN + R. (11) 

Since in a stable system in steady state the average number of messages received in a waiting time is equal to the average 
number of arrivals in a waiting time (a variation of Little's law) we have N = XW = X(Wd + W s ). Substituting this in (fTTI) 
we obtain 

As 2 

W s = sX{W d + W s ) + —. 

3 Note that (\\U\\ — r*)+ represents max(0, \\U\\ — r*) and U is a uniformly distributed random variable over the network region 1Z. 
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This implies 



W s = -?-W d + -^-^, (12) 
1-p 2(1 - p) 



Substituting (TO) and O in © yields ©. ■ 
In addition to the average waiting time of a classical M/G/l queue given in 0, the queueing delay also increases due to the 
collector's travel. Note that the E[(||Z7||-r*) + ] term can be further lower bounded by E[||U"||] —r*, where E[||f7||] = 0.383VA 

m- 

V. Collector Policies 

We derive upper bounds on delay via analyzing policies for the collector. The TSPN policy analyzed in Section IIII-B1I is 
stable for all loads p < 1 and has O(j^) delay scaling. Since the lower bound in Section [TV] also scales with the load as 

— the TSPN policy has optimal delay scaling. 

A. First Come First Serve (FCFS) Policy 

A straightforward policy is the FCFS policy where the messages are served in the order of their arrival times. A version of 
the FCFS policy, call FCFS', where the receiver has to return to the center of the network region (the median of the region 
for general network regions) after each message reception is shown to be optimal at light loads for the DTRP problem J6). 
This is because the center of the network region is the location that minimizes the expected distance to a uniformly distributed 
arrival. Since in our system we can do at least as good as the DTRP, FCFS' is optimal also for our system at light loads. 
Furthermore, the FCFS policy is not stable for all loads p < 1, namely, there exists a value p such that the system is unstable 
under FCFS policy for all p > p. This is because in the FCFS system the average travel component of the service time is 
fixed, which makes the average arrival rate greater than the average service rate as p — > 1. Therefore, it is better for a policy 
to serve more messages in the same "neighborhood" in order to reduce the amount of time spent on mobility. 

B. Partitioning Policy 

Next we propose a policy based on partitioning the network region into subregions and the collector performing a cyclic 
service of the subregions. This policy is an adaptation of the Partitioning policy of J6) to the case of a system with wireless 
transmission. We explicitly derive the delay expression for this policy and show that it scales with the load as O(j^) as in 
the TSPN policy. 

We divide the network region into (\/2r* xV2r*) squares as shown in Fig. [2] This choice ensures us that every location in 
the square is within the communication distance r* of the center of the square. The number of subregions in such a Partitioning 
is given bj^ n s = Aj (2(r*) 2 ). The partitioning in Fig. |2]represents the case of n s = 16 subregions. The collector services the 
subregions in a cyclic order as displayed in Fig. [2] by receiving the messages in each subregion from its center using an FCFS 
order. The messages within each subregion are served exhaustively, i.e., all the messages in a subregion are received before 
moving to the next subregion. The collector then receives the messages in the next subregion exhaustively using FCFS order 
and repeats this process. The distance traveled by the collector between each subregion is a constant equal to \[2r* . It is easy 
to verify that the Partitioning policy behaves as a multiuser M/G/l system with reservations (see |5j p. 198]) where the n s 

4 If y/ri7 is not even or if the subregions do not fit the network region for a particular choice of r*, then one can partition the region using the largest 
reception distance r* < r* such that these conditions are satisfied. 
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Fig. 2. The partitioning of the network region into square subregions of side \/2r* . The circle with radius r* represents the communication range and the 
dashed lines represent the collector's path. 

subregions correspond to users and the travel time between the subregions corresponds to the reservation interval. Using the 
delay expression for multiuser M/G/l queue with reservations in p. 200] we obtain, 



T, 



As 2 



part 



2(1 -p) 2v(l-p) 



V2r 



s, 



(13) 



where p = As is the system load. Combining this result with ([8]l and noting that the above expression is finite for all loads 
p < 1, we have established the following observation. 

Observation 1: The time average delay in the system scales as Q(jz^) with the load p and the Partitioning policy is stable 
for all p < 1. 

Despite the travel component of the service time, we can achieve 0( yz— ) delay as in classical queuing systems (e.q., G/G/l 
queue). This is the fundamental difference between this system and the corresponding system where wireless transmission is 
not used, as in the latter system the delay scaling with load is 9( (il p yi ) @- This difference can be explained intuitively as 
follows. Denote by N the average number of departures in a waiting time. It is easy to see from the P-K formula that in a 
classical M/G/l queue, N scales with the load as We argue that this scaling for N is preserved in our system but not 

in J6). The W s expression as a function of Wd in ( fT2b implies that for any given policy with its corresponding Wd, N can be 
lower bounded by yz- 4 . For the system in (6), the minimum per-message distance the collector moves in the high load regime 

scales as f2(^=) [6] (intuitively, the nearest neighbor distance among N uniformly distributed points on a square region of 

area A scales as ^=)- Therefore, for this system we have Wd ~ NQ(^=) rj f2(V NA) which gives N w il( j^z^p )■ Namely, 
Wd increases with the load and this results in an extra 1/(1 — p) scaling in delay in addition to the 1/(1 — p) factor of classical 
M/G/l queues. However, with the wireless reception capability, the collector does not need to move for messages that are 
inside a disk of radius r* around it. Since a finite (constant) number of such disks cover the network region, Wd can be upper 
bounded by a constant independent of the system load (for the Partitioning policy an easy upper bound on Wd is the length 
of one cyclic tour around the network). Therefore, in our system TV scales as 1/(1 — p) as in classical queues. 

In ifTSl we analyzed the case where the messages were transmitted to the collector using a random access scheme (i.e., 
with probability p in each time slot) and obtained S7( (jz^p ) delay scaling as in the system without wireless transmission. The 
reason for this is that in order to have successful transmissions under the random access interference of neighboring nodes, 
the reception distance should be of the same order as the nearest neighbor distances ITOl . lfT4l . 
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Fig. 3. Delay lower bound vs. network load using different SNR values for A = 200, 
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Fig. 4. Delay in the Partitioning policy vs the delay lower bound for SNR C = 17dB (r* = 2.2), 
(A = 800, v = 1, s = 2). Case-2: Comparable travel and reception times (A = 60, v = 10, s = 2). 



2 and o = 4. Case-1: Dominant travel time 



C. Numerical Results-Single Collector 

Here we present numerical results corresponding to the analysis in the previous sections. We lower bound the delay expression 
in © using E[(||(7|| - r*)+] > E[||i7||] - r* (where [||J7||] = 0.383-/A is the expected distance of a uniform arrival to the 
center of square region of area A (6]). Fig. [3] shows the delay lower bound as a function of the network load for different levels 
of channel qualitjO- As the channel quality increases, the message delay decreases as expected. For heavy loads, the delay in 
the system is significantly less than the delay in the corresponding system without wireless transmission in [6|, demonstrating 
the difference in the delay scaling between the two systems. For light loads and more noisy communication channels, the delay 
performance of the wireless network tends to the delay performance of (6). 

Fig. [4] compares the delay in the Partitioning Policy to the delay lower bound for two different cases. When the travel 
time dominates the reception time, the delay in the Partitioning policy is about 10.6 times the delay lower bound. For a more 
balanced case, i.e., when the reception time is comparable to the travel time, the delay ratio drops to 2.4. 

VI. Multiple collectors 

In this section we extend our analysis to a wireless network with multiple identical collectors. An arriving message is 
transmitted when one of the m collectors comes within the reception distance of the message location and grants access for 
the message's transmission. Therefore, at a given time there can be at most m transmissions in the network. We consider 

5 For the delay plot of the no-communication system, the point that is not smooth arises since the plot is the maximum of two delay lower bounds proposed 
in GO. 
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policies that partition the network region into m subregions. Each collector is assigned to one of the subregions and is allowed 
to operate only in its own subregion. We call these policies the class of network partitioning policies. In such a case, there 
is no interference from nodes within the subregion where the transmission is taking place. The only source of interference 
can be due to transmissions in other subregions. However, we assume that signaling schemes used in different subregions are 
orthogonal to each other so that there is no destructive interference from other subregions. This can be achieved for example 
by having different frequency bands for transmissions in different subregions. Furthermore, for networks deployed on a large 
area, even if orthogonal signaling is not utilized, the interference between subregions can be negligible due to signal attenuation 
with distance. Hence, we assume that simultaneous transmissions do not interfere with each other. Note that this is consistent 
for the purposes of a lower bound on delay since the interference from neighboring nodes can only increase the message 
reception time. Under these assumptions, utilizing the SNR criteria for successful message reception as in Section IIV1 the 
optimal reception distance is given by r* = (SNRc/^) 1 /" and each reception takes time s. 

A. Stability 

Here we show that p = Xs/m < 1 is a necessary and sufficient condition for stability of the system. 

1) Necessary Condition for Stability: A necessary condition for stability of the multi-collector system is given by p = 
Xs/m < 1. We prove this by showing that the system stochastically dominates the corresponding system with zero travel times 
(i.e., an M/D/m queue, a queue with Poisson arrivals, constant service time and m servers) similar to Section HIT1 

Theorem 4: A necessary condition for the stability of any policy is p = Xs/m < 1. Furthermore, the optimal steady state 
time average delay is lower bounded by 

T *> ^ H^if! + S (14) 

where p = Xs/m is the system load. 

The proof is similar to the proof of Theorem [T] and is given in Appendix- A. It makes use of the fact that the steady state time 
average delay in the system is at least as big as the delay in the equivalent system in which travel times are considered to be 
zero (i.e., v = oo). 

2) Sufficient Condition for Stability: We argue that p < 1 is also sufficient for stability. This can be seen by dividing the 
network region into m equal subregions, and performing a single-collector TSPN policy in each subregion. Since the arrival 
process is Poisson, each subregion receives an independent Poisson arrival process of intensity A / m. Furthermore, each collector 
performs a TSPN policy independently of the other collectors. Therefore, using the stability result of the single-collector TSPN 
policy, the systems in each subregion are stable if p < 1. We state this fact in the following theorem: 

Theorem 5: The system is stable under the multi-collector TSPN policy for all loads p = Xs /m < 1. 
Note that a similar delay analysis to the single-collector TSPN case shows that the multi-collector TSPN policy has O(yz^) 
delay scaling with the load p. 

B. Delay Lower Bound 

In addition to the delay lower bound given in (TPfl) . we provide another lower bound for the optimal delay and take 
their simple average. The following lemma states the second lower bound on delay. It is based on the convexity argument that 
when travel component of the waiting time is lower bounded by a constant, the equal area partitioning of the network region 
minimizes the resulting delay expression out of all area partitionings. 
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Theorem 6: For the class of network partitioning policies, the optimal steady state time average delay is lower bounded 

by 

T m > h s, (15) 

1 — p v 

where p = As/to is f/ie system load. 

Proof: Here we use an approach similar to the proof of Theorem[3] We divide the average delay T into three components: 

T = W d + W s + s. (16) 

The below lemma provides a bound for Wd, the average message waiting time due to the collectors's travel, using a result in 
lfl6l for the m-median problem. 
Lemma 2: 



W d > 1 3V?n " '-. (17) 

v 

Proof: Let SI be any set of points in di with |fi| = m. Let U be a uniformly distributed location in 5? independent of f2 
and define Z* = min„ S Q || U — v ||. Let the random variable Y be the distance from the center of a disk of area A/m to a 
uniformly distributed point within the disk. Then it is shown in ifTBI that 

E[f(Z*)] > E[f(Y)] (18) 

for any nondecreasing function /(.). Using this result we obtain E[max(0, Z* —r*)] > E[max(0, Y — r*)). Note that Wd can 
be lower bounded by the expected distance of a uniform arrival to the closest collector at the time of arrival less r* . Because 
the travel distance is nonnegative, we have 

W d > E[max(0,y-r*)]/« > max(0,E[F] -r*)/v, 



where the second bound is due to Jensen's inequality. Substituting E[Y] = fy^; into the above expression completes the 

proof. ■ 
Intuitively the best a priori placement of m points in 5? in order to minimize the distance of a uniformly distributed point 
in the region to the closest of these points is to cover the region with m disjoint disks of area A/m and place the points at 
the centers of the disks. Such a partitioning of the region is not possible, however, using this idea we can lower bound the 
expected distance as in dT8"| >. 

We now derive a lower bound on W s . Let R 1 , R 2 , R m be the network partitioning with areas A , A 2 , A m respectively 
(2~2jLi A? = A). Consider the message receptions in steady state that are received by collector j eventually. Let A J be the 
fraction of the arrival rate served by collector j. Due to the uniform distribution of the message locations we have 

X j _ AJ 

T ~ T ' 

Let be the average number of message receptions for which the messages that are served by collector j waits in steady 
state. Similarly let W% and W° d be the average waiting times for messages served by collector j due to the time spent on 
message receptions and collector j's travel respectively. Using ( fTTT i and lower bounding the residual time by zero we have 

Wl > sN j . 
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Using Little's law (N J = X 1 (Wj> + Wi)) similar to the derivation of (12\ we have 



wi > AJ f. wi. (i9) 

s - 1 - A^s d 



The fraction of messages served by collector j is A> /A. Therefore, we can write W s as 

m A j 

3=1 

\ W/3 



> > -r-^ — r— Wi- (20) 



3 = 1 



For a given region B? with area A\ Wi is lower bounded by (similar to the derivation of ( fTOl i) the distance of a uniform 
arrival to the median of the region less r*. 



^ E[max(0, \\U - v\\ - r*)] 



v 



> max(0,E[||t/-H|]-r^ 
v 

where v is the median of R? and \\U — u\\ is the distance of U, a uniformly distributed location inside RP, to i/. The inequality 
in (|2TT > is due to Jensen's inequality for convex functions. A disk shaped region yields the minimum expected distance of a 
uniform arrival to the median of the region. Using this we further lower bound Wa by noting that for a disk shaped region of 

area Aj, E[||C7 — is just the expected distance of a uniform arrival to the center of the disk given by — ■ Hence 



W j > max(0,f^ff-r*) = ma x(0, Cl - r*) ^ 



where Cl = ^h= = 0.376. Letting f(A?) = , which is a convex and increasing function of A\ we rewrite d20l l as 

v 1r 1— \ — s 

m f(A^) 

W s > V ^-H-^max(0, Ci Va1 - r*). (23) 

3 = 1 



Next we will show that the function /(A-')A : 'max(0, CiV A? — r*) is a convex function of A 3 via the two lemmas below. 

Lemma 3: Let /(.) and <?(.) fee fwo convex and increasing functions (possibly nonlinear) defined on [0, A]. The function 
h(.) = f.g(.) is also convex and increasing on its domain [0, A]x[0, A]. 

Proof: See Appendix-B. ■ 
Lemma 4: h(x) = x max(0, C\ \fx — Ci) is a convex and increasing function of x. 

Proof: See Appendix-C. ■ 
Letting g{A>) = f(A J )A J max(0, c\\f~AP — r*), we have from the lemmas [3] and H] that the function g(A 3 ) is convex. Now 
rewriting ( 1231 we have 

3=1 
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Using the convexity of the function g(A^) we have 




) 



m. 




|f max(0,ciVI7 



m — r*) 



As 

m 



/0 




(24) 



1 



The above analysis essentially implies that the W s expression in (l23l l is minimized by the equitable partitioning of the network 



Finally, taking the simple average of $15[ and ( TPfl i we arrive at the following theorem. 

Theorem 7: For the class of network partitioning policies, the optimal steady state time average delay T* n is lower bounded 



where p = Xs/m is the system load. 

Theorem [7] is valid for the class of network partitioning policies. For the system without wireless transmission, it has been 
shown that partitioning the region into m equal size disjoint subregions (one for each collector) preserves optimality in the 
high load limit Q, BP . We conjecture that this optimality is also preserved in our system. 

C. Multiple Collector Policies 

Note that it is shown in Q that a generalization of the FCFS policy, namely, creating m Voronoi regions with centers of the 
regions given by the ?7i-median location^] of the network region and having each receiver perform the single-receiver FCFS 
service in each region has optimal delay at light loads. This also holds in our system via a similar argument to Section IVl 
Moreover, this policy is not stable as p — > 1 due to the same reason as in Section [V] 

1) Generalized Partitioning Policy: Next we propose a policy based on dividing the network region into m equal size 
subregions. Each collector is assigned to one of the subregions and is responsible for receiving messages that arrive into its 
own subregion using the single collector partitioning policy analyzed in Section [V-BI Namely, first the network region is divided 
into subregions of area A/m and then each subregion is divided into \[2r* x\/2r* square j^. The number of \/2r* xy/2r* 
squares in each subregion is given by n s = ■ ^8- represents such a partitioning for the case of four collectors in the 

network with n s = 16 squares in each subregion. Since each subregion behaves identically, the average delay of this policy 
is the average delay of the single collector Partitioning policy applied to a subregion with arrival rate A/m, area A/m, and 

6 The set of m-median locations for a region is the set of the best m a priori locations in the region that minimizes the expected distance to a uniform 
arrival. 

7 Note that the number of collectors must be a square number in order to divide the network region into subregions of exactly equal shape and size. 



region. Finally combining ( fT6l ), (JT7J and (l24t we obtain (115) . 



by 




(25) 
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Fig. 5. The partitioning of the network region into square subregions of side \/2r* for the case of multiple collectors in the network. The circle with radius 
r* represents the communication range and the dashed lines represent one of the collector's path. 
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Fig. 6. Delay lower bound vs. network load for m=2 collectors, SNR C = 30dB (r* = 4.7), A = 400, ft = 2, a = 4, v = 1 and s = 1. 



A/m 

n s — 2 (r*)2 : 



1 part 



2m(r*) 2 



2m(l - p) 2u(l - p) 



(26) 



where p = As/to is the system load. This result, when combined with ( fl4l i. establishes that for the case of multiple collectors 
in the system the delay scaling with the load is 9(— ^). This is again a fundamental improvement compared to the 6( nj^ys ) 
delay scaling in the system without wireless transmission and with multiple collectors in J7J. 

2) Numerical Results: We compare the delay lower bound in d25l l to the delay lower bound in the corresponding system 
without wireless transmission in Q for the case of two collectors in Fig. [6] The delay in the two-collector system is significantly 
below the delay in the system without wireless transmission and this difference is more pronounced for high loads. Fig. [7] 
displays the delay lower bound in (l25"l l and the delay of the Partitioning policy in d26l l as functions of the network load p. The 
delay of the Partitioning policy is about 7 times the delay lower bound. 

VII. CONCLUSION 

In this paper we considered the use of dynamic vehicle routing in order to improve the delay performance of wireless 
networks where messages arriving randomly in time and space are gathered by a mobile collector. We characterized the 
stability region of this system to be all system loads p < 1 and derived fundamental lower bounds on time average expected 
delay. We derived upper bounds on delay by analyzing policies and extended our results to the case of multiple collectors in 
the system. Our results show that combining controlled mobility and wireless transmission results in 0(j3^) delay scaling 
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Fig. 7. Delay of the Partitioning policy vs the delay lower bound for m = 4 collectors, SNR C = 20dB (r* = 2.6), A = 500, /3 = 2, a = 4, v = 1 and 
s = 2. 

with load p. This is the fundamental difference between our system and the system without wireless transmission (DTRP) 
analyzed in J6) and Q where the delay scaling with the load is 6( ( 1 _ 1 p )2 )■ 

This work is a first attempt towards utilizing a combination of controlled mobility and wireless transmission for data collection 
in stochastic and dynamic wireless networks. Therefore, there are many related open problems. In this paper we have utilized 
a simple wireless communication model based on a communication range. In the future we intend to study more advanced 
wireless communication models such as modeling the transmission rate as a function of the transmission distance. Finally, 
extending our results for a general message location distribution in the network is a subject of future research. 

Appendix A-Proof of Theorem[5] 
The proof is similar to the proof of Theorem Q] First consider the following lemma. 

Lemma 5: The steady state time average delay in the system is at least as big as the delay in the equivalent system in which 
travel times are considered to be zero (i.e., v = oo). 

Proof: The proof is similar to the proof of Theorem[5] Consider the summation of per-message reception and travel times, 
s and di, as the total service requirement of a message in each system. Since di is zero for all i in the infinite velocity system 
and since the reception times are constant equal to s for both systems, the total service requirement of each message in our 
system is deterministic ally greater than that of the same message in the infinite velocity system. Let Di,D2, ■ •■,-D n be the 
departure instants of the first, second and similarly the n th message in the system. Similarly let D 1 , D 2 , D n be the departure 
instants of the first, second and similarly the n th message in the infinite velocity system. Similarly let Ai,A2, A n be the 
arrival times of the first second and similarly the n th message in both systems. We will use complete induction to prove that 
Di > D i for all i. Since the service requirement of each message is less in the infinite velocity system, we have D\ > D v 
Assume we have Di > D i for all i < n. We need to show that D n+ i > D n+1 in order to complete the complete induction. 
We have 

A n +i < D n+ i - s, (27) 
hence the n + 1 th message is available at time D n+ i — s. We also have 

"n+l— m — ^n+l-m < -Dn+1 — S. 

The first inequality is due to the complete induction hypothesis and the second inequality is due the fact that the m th last 
departure before the n + 1 th departure has to occur before the time D n+ i — s. Hence there is at least one collector available in 
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the infinite velocity system before the time D n+ i — s. Combining this with (|27l proves the complete induction. Now let D(t) 
and D (t) be the total number of departures by time t in our system and the infinite velocity system respectively. Similarly let 
N(t) and N (t) be the total number of messages in the two systems at time t. Finally let A(t) be the total number of arrivals 
by time t in both systems. We have N(t) = A(t) — D(t) and ./V (t) = A(t) — D (i). From the above induction we have 
D(t) < D (t) and therefore N(t) > N (t). Since this is true at all times, we have that the time average number of customers 
in the system is greater than that in the infinite velocity system. Finally using Little's law proves the lemma. ■ 
When the travel time is considered to be zero, the system becomes an M/D/m queue (a queue with Poisson arrivals, constant 
service time and m servers). Therefore we can bound using bounds for general G/G/m systems. In particular, the waiting 
time Wa/G/m m a G/G/m queue with service time s is bounded below by 0181 p. 48] 

A 777 — 1 

W G/G/m >W — , (28) 

where W is the waiting time in a single server system with the same arrivals as in the G/G/m queue and service time s/ra. 
Since in our case the infinite velocity system behaves as an M/D/m system, W has an exact expression given by the P-K 
formula: W = As 2 /(2m 2 (l — p)) where p = Xs/m. Substituting this in (|28| | and using Lemma|5]we have ( fT4l . 

Appendix B-Proof of Lemma[3] 
Clearly h(.) = f.g(.) is increasing. Let x and y be two points in the domain of h and let a G (0, 1) be a real number. 

h(ax + (1 — a)y) = f(ax + (1 — a)y)g(ax + (1 — a)y) 
< (af(x) + (1 - a)f{y)). 

.(ag(x) + (1 - a)g(y)) 
= a 2 f(x)g(x) + (1 - a) 2 f(y)g(y) 

+a(l-a)f(x)g(y)+a(l-a)f(y)g(x)), 

where the inequality is due to the convexity of / and g. We add and subtract af(x)g(x) and after some algebra obtain 

h(ax + (1 — a)y) < ah{x) + (1 — a)h(y) 

+a(l-a)(f(x)- f(y))(g(y)-g(x)) 
< ah(x) + (1 - a)h(y), 

where the last inequality is due to the fact that / and g are increasing functions. 



17 



Appendix C-Proof of Lemma[4] 

It is clear that h(x) is an increasing function of x. Let x and y be two points in the domain of h and let a <G (0, 1) be a 
real number. 

h(ax +(1 — a)y) = 

= (ax + (1 — a)y)max(0, c\ y/ ax + (1 — a)y — C2) 
= max(0, ci(ax + (1 — a)y)? — c 2 (ax + (1 — a)y)) 

< max(0, c\(ax^ + (1 — a)y^) — C2(cke + (1 — a)y)) 
= max(0,aa;(ci^ - c 2 ) + (1 - a)y(c ly /y - c 2 )) 

< max(0,aa;(ci^-C2))+max(0, (l-a)y(ci^y-c 2 )) 
= ah(x) + (1 — a)h(y), 

where the first inequality is due to the convexity of the function x 2. This shows that h(x) is a convex and increasing function. 
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